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ABSTRACT 

Direct detection of thermal emission from nearby hot Jupiters has greatly 
advanced our knowledge of extrasolar planets in recent years. Since hot Jupiter 
systems can be regarded as analogs of high contrast binaries, ground-based in- 
frared long baseline interferometers have the potential to resolve them and detect 
their thermal emission with precision closure phase - a method that is immune 
to the systematic errors induced by the Earth's atmosphere. In this work, we 
present closure phase studies toward direct detection of nearby hot Jupiters us- 
ing the CHARA interferometer array outfitted with the MIRC instrument. We 
carry out closure phase simulations and conduct a large number of observations 
for the best candidate v And. Our experiments suggest the method is feasible 
with highly stable and precise closure phases. However, we also find much larger 
systematic errors than expected in the observations, most likely caused by disper- 
sion across different wavelengths. We find that using higher spectral resolution 
modes (e.g., R=150) can significantly reduce the systematics. By combining all 
calibrators in an observing run together, we are able to roughly re-calibrate the 
lower spectral resolution data, allowing us to obtain upper limits of the star- 
planet contrast ratios of v And b across the H band. The data also allow us to 
get a refined stellar radius of 1.625±0.011 R©. Our best upper limit corresponds 
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to a contrast ratio of 2.1x10^:1 with 90% confidence level at 1.52/im , suggesting 
that we are starting to have the capability of constraining atmospheric models 
of hot Jupiters with interferometry. With recent and upcoming improvements 
of CHARA/MIRC, the prospect of detecting emission from hot Jupiters with 
closure phases is promising. 

Subject headings: infrared: stars - planets: hot Jupiters - planets: v And b - 
techniques: interferometry - facility: CHARA 



Introduction 



The discovery of a planet ar ound a nearby star 51 P eg in 1995 opened a window into new 
worlds outside the solar system (iMayor &: Quelozlll995h . Since then, more than 500 so-called 
exoplanet^ have been discovered, revolutionizing our knowledge of their nature and origin. 
Among those discoveries, about 24 planets had their thermal emission directly detected by 
photometric and/or spectroscopic measurements from space or grou r id, most of which are 



(e.g. 



Beerer et al. 2010: Machalek et al. 2010 



known as "hot Jupit ers" or "hot Neptunes' 
Nymeyer et al.ll2010l . etc.). So far, direct detection of thermal emission and characterization 
of planetary atmospheres is generally only possible for these hot planets as they are very 
close to their host stars (<0.1AU), and are thus heated enough to have temperatures above 
lOOOK, providing as high as 10~^ of their host stars' flux in the infrared (e.g., J, if, and 
mid-infrared) . 



The atmospheres of hot Jupiters have many interesting properties. For instance, trans- 
missio n spectra at primary eclipses have show n the presence of sodium, water and methane 



(e.g., iRedfield et al.ll2008l : ISwain et al.ll2009l ). while thermal emission at se condary eclipses 
has shown the presence of water , carbon dioxide, and carbon monoxide (e.g.. lKnutsonll2007al : 
Barmanll2008l : ISwain et al.ll2008l ). Due to their close-in orbits, hot Jupiters are tidally locked 
to their host stars, leading to a constant day side that experie nces intense stellar irradi- 
ation and a cold night side that remains in perpetual shadow (IGuillot et al.lll996l ). The 
temperature difference between the day and night sid es thus induces atmospheric circula- 
tion and strong zonal winds to redistribute the heat (iKnutsonI l2007bl ) . Studies have also 
suggested the existence of a hot stratosphere on the dayside of some hot Jupiters (e.g., HD 
209458b and HD 149026b), which inverts the te mperature profile of higher atmospher i c lay- 



ers, and flips water absorptions into emissions (iKnutson et al.l 12008 



Burrows et al. 2008 



data from The Extrasolar Planets Encyclopaedia: 



http://exoplanet.eu/catalog.php 
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Madhusudhan fc Seagerl 120101) . However, some other planets such as HP 1 89733b seem to 



lack such thermal inversion ( ICharbonneau et al.l l2008l : iGrillmair et al.l |2008| ) , implying fun- 



damental atn aospheric differences between these planets and leading to sub-classification of 
hot Jupiters (IBurrows et al.ll2008l ). Studies of hot Jupiters' atmospheres will not only reveal 
their composition, structure, and dynamics, but will also shed light on our understanding 
of the planet formation processes. More importantly, characterizing hot Jupiters allows us 
to pave the path toward characterizations of Super-Earth planets and eventually. Earth-like 
planets. 

Among the ~ 22 detected hot Jupiters, however, a majority of them are transiting 
planets whose orbits are aligned with the line of sight from the Earth, and only 2 of them are 
non-tr ansiting planets (i.e., v And b and HD 179949b. [Harrington et al.l ( 120061 ) : ICowan et al. 



( 120071 )). The lack of studies of non-transiting hot Jupiters leaves a great opportunity for 
long baseline optical/infrared interferometry, in that long baseline interferometers can see hot 
Jupiter systems as extremely high-contrast binaries, and thus can directly determine their 
orbital elements and provide accurate mass estimates. Interferometric measurements can also 
provide absolute planet/star flux ratios (or star/planet contrasts) of non-transiting planets 
which cannot be disentangled by combined-light techniques (i.e., through transits). In fact, 
interferometric measurements also provide an independent way of characterizing transiting 
planets in addition to the combined light techniq ue. Since the bulk of energy from hot 
Jupiters emerges from the near-IR between 1-3 /im ( IBurrows et al.ll2008l ). interferometry can 
provide a better understanding of their spectra and global energy budget with measurements 
at near-IR bands. 

However, detecting the weak emission from a hot Jupiter from the ground is a difficult 
task. To date, only a few hot Jupiters' thermal emission has been detected from the ground, 
while the rest were all achieved from space. To reach this goal, we require very stable and 
high precision measurements, and most importantly, require a method that can eliminate 
or calibrate the effects of Earth's atmospheric turbulence. One possible approach is to use 
high precision and high resolution closure phase measurements obtained with ground-based 
interferometers. Closure phase is measured by combining the phases of three baselines in 
a closed triangle. It is immune to any phase shifts induced by the atmospheric turbulence, 
including the differential chromatic dispersion that affects the differential phase. The major 
biases or systematic errors of closure phase come from non-closed triangles introduced in 
the measurement process, and in principle, can be precisely calibrated. Therefore, it is a 
good observable for stable and preci s e measurern ents. More descriptions of the closure phase 
technique can be found in iMonnierl ( l2003l . 120071 ). 



Studies have already been carried out to explore the possibilities of using closure phases 
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for exopl anet detection (e.g ., IJoergens &: QuirrenbachI |2004J : IZhao et al.ll2008al . etc.). Par- 



ticularly, IChelli et al.l ( 120091 ) studied in detail the characteristics of closure phases and the 
corresponding SNR when the primary star of a binary is getting r esolved and approac hing 
the visibility null (i.e, "Phase Closure Nulling"). Using this method. iDuvert et al. fcoioh de- 
tected th e faint close compa nion of HD 59717 at a contrast of ~ 100 : 1 using VLTI/AMBER. 



Recently, lAbsil et al.l ( l2010l ) also excluded the presence of a brown-dwarf companion in the 
innermost region of the /5 Pic planetary system at a upper contrast limit of ~ 200 : 1 using 
closure phases obtained by VLTI/AMBER. 

In this paper, we report our closure phase studies toward the ultimate goal of directly 
detecting emission from hot Jupiters using CHARA/MIRC. The paper is organized as follows. 
We first briefly introduce our candidate v And b and our observations in §21 In §3]we simulate 
the closure phase signals and the required precision for the candidate. We then discuss the 
calibration issues we encountered in test observations and present our solutions. Based on 
our calibrations, we present a preliminary upper limit for v And b in ^ Finally, we conclude 
our studies and give future prospects in ^ 



2. Candidates and observations 



2.1. Candidates 



Among all the known non-transiting hot Jupiters, several of them are within the sensi- 
tivity limit of CHARA/MIRC. v And b is currently the most favorable candidate due to its 
relatively high temperature and the high brightness of its host star. Therefore, in this study 
we focus only on this best candidate. 



V And is an F8V star located 13.5 pc away from the Sun. iButler et al.l (Il997| ) first 
discovered its hot Jupiter v And b in 199 7, which has a perio d of 4.6 days and is orbiting 
at 0.06 AU. The follow-up observations of iButler et al.l (Il999[ ) found two more companions 
in the system, v And c and v And d. v And c is orbiting at 0.83 AU from the host star 
with a period of 24 1 days, while v And d is orbiting at 2.5 AU with a period of 1267 days 
( IButler et al.lll999l ). Most recently, ICuriel et al.l (120 111 ) found a fourth planet v And e with a 
period of 3848.9 days at 5.25 AU. The syste m is non-coplanar. Th e inclination of the closest 
planet v And b is found to be likely > 28° (ICrossfield et al.ll2010l ). while t he second and the 



third planets {v And c & d) have a relative inclination of 15°-20°. In 2006, iHarrington et al. 



( I2OO6I ) directly detected thermal emission from the hot Jupiter v And b using Spitzer MIPS 
at 24/xm, in which they detected the relative day-night flux variations of the planet over five 
epochs of the whole 4.6-day orbital period, and provided a lower limit to the planet/star 
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flux ratio. Later, ICrossfield et al.l ( 120101 ) further refined the flux variation curve with more 
MIPS 24/im data, and found an unusually large phase shift of the flux maximum (~ 80"). 
Atmospheric models have also been applied to interpret these Spitzer data. However, as 



Burrows et al.l (120081 ) pointed out, due to the lack of absolute flux level and information in 



other wavelengths, there are too many degrees of freedom to draw strong conclusions about 
the planetary and atmospheric properties. Thus, detection of its absolute planet/star flux 
ratios and at other wavelengths such as near-IR are required to break model degeneracies. 
Furthermore, high si gnal-to-noise detecti ons will even allow us to obtain the absolute phase 
curve of the planet (IBarman et al.ll2005l ). providing constraints to the circulation and heat 
redistribution patterns of its atmosphere. 



2.2. Observations 



Closure phase measurements require at least three telescopes to be combined in a closed 
triangle. Currently, the MIRC jMonnier et al.l booj ) and CLIMB ( Isturmann et al.l I2OIC 
instruments at the CHARA array, and the AMBER JPetrov et al.l I2OO7I ) and PIONIER 
instruments at the VLTI have this capability in the near-IR. In this study, we employ 
CHARA/MIRC for our measurements. 



The CHARA array, located on Mt. Wilson, consists of six 1-meter telescopes (jten Brummelaar et al. 



2OO5I ). The array is arranged in a Y-shaped configuration to provide good position angle 
coverage. The six telescopes form 15 baselines ranging from 34m to 331m, making CHARA 
the longest-baseline optical/IR interferometer array in the world and providing a resolution 
of ~0.5 mas in the H band. The Michigan Infra-Red Combiner (MIRC) is an imaging com- 
biner that currently combines 4 CHARA telescopes, providing 6 visibilities, 4 closure phases 
and 4 triple amplitudes simultaneously. MIRC works at both H and K bands, and has three 
spectral modes: R=40 (with prism), 150, & 500 (with grism). The lowest resolution mode 
(R=40) disperses light into 8 spectral channels on the detector, while th e R=150 and R=50(3 
modes disperse light into 24 channels and 80 channels respectively (see lMonnier et al.ll2004j . 



20061 . for details). The compact design of MIRC allows for stable calibration and precise 



closure phase measurements, thus it is best suited for the purpose of this study. 

We conducted observations of v And on 23 nights from 2006 to 201 with a total 



integ r ation time of 8.2 hours, following the standard observing procedures (jMonnier et al. 



2007 



Zhao et al.l |2009[ ) . The observation log is listed in Tabled! Several combinations 



of four CHARA telescopes are used in the observations, while in most cases we adopt the 



^http://www-laog.obs.ujf-grenoble.fr/twiki/bin/view/Ipag/Projets/Pionier/WebHonie 
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combination S1-E1-W1-W2 and S2-E2-W1-W2 for good baseline coverage. The observations 
were mostly conducted in the H band (A=1.5 - 1.8 yum) with the lowest resolution mode of 
R=40, except for two nights in 2010 when we employed the higher resolution mode of R=150. 
Each observation of the target was bracketed with calibrators for visibility and closure phase 
calibration. For the purpose of bias subtraction and flux calibration, each set of fringe data is 
bracketed with measurements of background (i.e., data taken with all beams closed), shutter 
sequences (i.e., data taken with only one beam open at a time to estimate the amount of light 
coming f rom each bearn), and foreground (i.e., data taken with all beams open but without 
fringes) (jPedretti et al.ll2009[ ). Each object is observed for multiple sets. Du ring the period 



of tak ing fringe data, a group-delay fringe tracker is used to track the fringes ( iThureau et al. 
20061 ). In order to track the flux coupled into each beam in "real time" to improve the 



visibility measurements, we used spinning choppers to temporally modulate the light going 
into each fiber simultaneo usly with fringe m easurements. In 2009, photometric channels were 
commissioned for MIRC (IChe et al.l l2010l ). and the choppers were made obsolete because 
of the much better real-time flux calibration provided by the photometric cha nnels. The 
data reduction process also follows the pipeline outlined in (IMonnier et al.l 120071 ). where the 
closure phases are extracted from the phase term of the complex triple amplitudes after the 
subtraction of background, foreground, and correction of fiber coupling efficiencies. 



3. Closure phase simulations and calibration studies 
3.1. Simulations 



Because exoplanet-host stars and their hot Jupiters are similar to high contrast close 
binaries, we simulate the closure phase signals using binary models for v And b. We choose 
to use the longest telescope triangle of the CHARA arra y (i.e., Sl-El-Wl) in our study to 
obtain the highest closure phase levels (IChelli et al.ll2009l ). The latest orbital properties of v 
And b and the size of its host st ar are listed in Table El The Infrared planet/star flux ratios 
are adopted from the models of ISudarsky et al.l (120031 ). 



Figure [T] shows the closure phase simulations for v And b at four wavelength channels. 
Since the value of the ascending node [Vt) is unknown for v And b, we assume 4 different 
values in our simulations: VL = 0°, 45", 90", 135". The inclination is fixed to the most probable 
values assuming that the planet is coplanar with the star's rotation (see Table |2]). Figure [1] 
shows that the closure phase signal varies rapidly as the baseline projection varies during the 
Earth's rotation. Shorter wavelength channels generally have higher signal amplitudes and 
peak at ~ ±0.08°, due to the fact that the host star is more resolved at those wavelengths. 
The overall closure phase level of v And b is insensitive to the values of Q in the simulation 
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because of CHARA's even baseline layout (i.e., the Y-shape), which covers all position angles 
more or less equally well. We note here that because the visibility of the star v And only 
goes to null at certain wavelength channels, the closure phase signal in the simulation is very 
sensitive to the value of adopted stellar diameter. Figure [T] also indicates that in order to 
detect the closure phase signal from v And b, the required 1-a precision has to be better 
than roughly ±0.05°. 



3.2. Calibration studies 

To compare the precision of our measurements with the simulations, we first examine 
the quahty and stability of our data. The left panels of Figure [2] show a good night of 
closure phase measurements for the middle wavelength channel of MIRC, obtained with 
the largest triangle of CHARA (Sl-El-Wl). The closure phases are stable over the 1.5 
hrs of observation and the error averages down roughly as a/ZV (the bottom left panel), 
suggesting the measurements are immune from systematic errors like changes in the seeing. 
The nominal measurement error is 0.3° when averaging the whole 1.5 hrs of observation 
together. The performance at a similar channel wit h a shorter t riangle of CHARA (E2-W1- 



W2) is demonstrated to be about 3 times better in IZhao et al.l (j2008bl ). 



Although these measurements are promising and stable, we have also encountered large 
unexpected systematic errors in other wavelength channels. As indicated in the right panels 
of Figure [21 the closure phases show a systematic drift and the associated errors cannot 
be averaged down as V^. In fact, the closure phases not only change with time, but also 
vary as a function of wavelength. More interestingly, the closure phase drifts are highly 
correlated for all calibrators within a whole observing run of many nights. Figure E] plots 
closure phases of six calibrators obtained in 6 nights in 2008 August (Algenib, 38 Tau, 7 
Lyr, 7 Tri, C, Peg, C Per) versus azimuth angles. There is a clear trend of closure phase 
change as a function of azimuth in the figure. Figure H] shows the similar correlation of 
closure phase with altitude. The closure phase drift reaches about 8 degrees between the top 
and the bottom panel. In addition, we can also see an obvious slope change centered at the 
middle wavelength channels in both Figures |3] & HI Although only 3 wavelength channels 
are shown here, the slope actually changes gradually from the first to the last channels of 
MIRC, and similar effects are also seen in other observing runs. Although the data shown 
in these figures were obtained in 6 nights, the correlations with altitude and azimuth are 
strong and consistent, suggesting that the major cause of the closure phase drifts may stem 
from the changing positions of the targets on the sky. 



Possible causes of the closure phase drifts may include: 1. polarization effects caused 
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by the non-identical beam trains of CHARA; and 2. extra dispersions in the delayhnes that 
are not compensated in vacuum, which can contaminate closure phases across wavelength 
channels and change with target position. To investigate the individual causes and find out 
the best solution, we carried out a series of experiments with: 1. using a linear polarizer 
to reduce the effect of polarization; 2. using a linear polarizer and a 40/im slit to reduce 
the effect of polarization and partially reduce dispersion; and 3. using a grism of R~150 to 
reduce dispersion only. We determine the slopes of the closure phase change as a function of 
only altitude for simplicity. Since the slopes of the closure phases also drift across wavelength 
channels (see Figure Hj), we use the magnitude of the slope change, i.e., (slope of channel 8 
— slope of channel 1), to characterize the closure phase changes. The magnitudes of slope 
change are averaged over four closure phase triangles of each observing run, and the errors 
are determined from the scatter of the four triangles. The results of the experiments are 
shown and compared in Table Ei 

Table [3] shows that the original observations have the largest magnitude of closure phase 
change. Observations with polarizer and polarizer + 40/im slit have slightly smaller magni- 
tudes of closure phase change. However, we have also seen some nights without polarizers 
have similar or even lower drifts than those with the polarizer, indicating the effect of using 
a polarizer is small and the major cause of the drifts may not be polarization effects. When 
using the grism of R=150, however, the correlated closure phase drifts and slope change 
become nearly zero. Although we only had a very small amount of data for the experiments, 
this preliminary result indicates that the extra air dispersion from the delaylines is most 
likely the major cause of the closure phase drifts, which contaminates closure phases with 
non-closed triangles from other wavelengths, especially at the edges of the bandpass. The 
dispersion effect also explains the slope change centered at the middle wavelength channels 
seen in Figures [3] & HI 

Since the closure phase changes shown in Figures E] & IH are highly correlated for all 
nights within an observing run with consistent system and optical settings. We can use all 
calibrators from a run to look for a solution to calibrate the drifts. We have experimented 
with several function forms to characterize the slopes, including linear function, linear sur- 
face, quadratic surface, etc. A linear surface fit can estimate the closure phase drifts well, 
while a quadratic surface fit works the best. Figure O shows an example of our best ap- 
proach of fitting the closure phases with quadratic surface functions of both target altitude 
and azimuth. The quadratic surface function characterizes the drifts very well for nearly all 
wavelength channels, as indicated by the improved values in the plot. It therefore can be 
employed as a new empirical model to calibrate our data within the same observing run. 

There is also a caveat, however, that this calibration scheme requires a wide span of 
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calibrator positions on the sky for a good coverage of altitude and azimuth to bracket the 
target, so that the quadratic fit can reliably estimate the closure phase change. However, 
most of our observations have a limited number of calibrators and calibrator visits, and 
thus do not have a wide range of position coverage. Therefore, we adopt the results from a 
linear surface fit for those cases instead in order to avoid unreliable or erroneous quadratic 
extrapolation, with a trade-off of slightly less accurate closure phase calibration. Further 
validations of these schemes and their robustness are required in a future work. 



4. Diameter and upper limits for v And b 
4.1. Refined diameter 



We first determine the angular diameter of v And with a uniform disk (UD) model and 
a limb-darkened disk (LD) model. The squared- v i sibilit ies {y"^) of the data are calibrated 
in the regular way as described in iMonnier et al.l (120071 ). UD diameters of the calibrators 
are listed in Tab le HI We apply a power law limb darkening, /(yu) = /i", for the LD model 
( iHestrofferl 119971 ) . Since the limb darkening coefficient a is similar for stars with the same 
spectral type, we adopt a fixed a and only vary the diameter in our L D fit . The value of a 
is interpolated and converted from the square root law coefficients of Ivan Hammd (I19931) . A 



value of a ~ 0.18 is estimated for v And, assuming [Fe/if]=0.15 fiButler et al.ll2006l ). We 



bootstraplfl the V"^ data from different nights to simulate the statistics. Each bootstrapped 
data set is then fit with a UD and LD model separately with minimization. A number 
of 150 bootstrap iterations are carried out and the median of the 150 best-fit diameters are 
taken as the global best-fit, while the l-cr error is determined from the scatter. Figure |6] shows 
the best-fit models of v And, together with the b inned data. Our final best-fit UD of v And is 
1.097±0.009 mas, consistent with the FLUOR ( iCoude du Foresto et al.ll2003l ) measurement 
of 1.0 98=b0.007 mas (Merand 2008, private communication), and the results of iBaines et al. 
J2008h . 1.09 1±0.009 mas. The best-fit LD size is 1.121±0.007 mas, also consistent with 



the result of iBaines et al.l (120081 ). and corresponding to a radius of 1.625±0.011 R© for the 



•^Bootstrapping is a technique that can provide robust simulations of the distribution of a data set. It is 
very useful for data sets with complicate d or unknown distributions, and is widely use d to derive estimates of 
standard errors and confidence intervals (jPress et al.lll992l : lEfron fc Tibshirani II 19931 ). Since we use multiple 
nights of data with various of baselines and different systematics (within each night) for joint solutions, the 
distributions of our parameters of interest are unknown. Thus, bootstrapping is a suitable technique for 
our data. Bootstrapping requires the assumption that the data are independent and identically distributed. 
This assumption holds for our data because we treat each night of data equally, each night's are independent 
and are acquired in the same way. 
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parallax of 74.12 ± 0.19 mas (Ivan Leeuweru 120071 ) (with error propagated). 



4.2. Upper limits 

Since we do not have enough SNR to detect the planet from a single night of observation, 
we can take the advantage of the well known orbital parameters of v And b and combine all 
observations in Table [1] together to increase the total SNR for higher precision. Due to the 
quickly varying closure phases caused by the Earth's rotation, as shown in Figured! we split 
the data into small chunks with an averaging time of less than lOmin to avoid smearing the 
signal. 

To calibrate the large closure phase systematics described in last section, we apply 
the new calibration scheme by fitting a quadratic surface function to all calibrators within 
each observing run, and subtracting the predicted zero closure phases from the raw values. 
For observing runs with inadequate calibrator coverages, we use the results from linear 
surface fit instead of quadratic fit. The total uncertainty of each data point is estimated 
by bootstrapping the data used in the quadratic or linear surface fit, and combining the 
additional uncertainty with the original values. Figure [7] compares the results before and 
after applying the new calibration scheme for the 23 nights of data. As we can see in Figure 
O the non-zero closure phases are calibrated out and the large scatters are reduced in the 
newly calibrated data. 

We then fit the new closure phase data with binary models to search for the closure 
phase signal from the planet. Because the closure phases are different for each channel, we 
determine planet/star flux ratios for each of the 8 wavelength channels simultaneously at 
the same orbital position, and search for the best position with a joint minimization. 
The orbital positions of the system are calculated by fixing the well known parameters from 
radial velocity observations, i.e., P, Tq, e, and u, and only varying the semi-major axis, Q, 
and inclination. The method and the corres ponding codes are validated with data of the 



well-known binary l Peg from iMonnierl (120071 ) . 



Although we have searched the parameter space extensively, we do not clearly detect v 
And b in our fits, suggesting that our current signal-to-noise ratio is still inadequate. We thus 
decided to report our results in terms of an upper limit to the planet/star flux ratio. To do 
this, we simulate the statistics of the best-fit planet/star flux ratios by bootstrapping different 
nights of data. This approach treats each night of data equally, ensuring the robustness of 
the bootstrap by preventing data from certain nights dominating the statistics. Because 
noise is dominant in the data, we use a fine grid search to ensure robust results. For each 
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bootstrapped data set, we searched an extensive range of semi-major axis, inclination, 
and the planet/star flux ratios at each of the 8 spectral channels. The set of parameters that 
yields the minimum is then chosen as the "best-fit" . A total of 150 bootstrap iterations 
are carried out, and the corresponding distributions of the best-fit planet /star flux ratios are 
shown in Figure [HI 

The dashed lines in Figure [8] indicate the upper limits of 90% confidence level, i.e., 
with 90% of chance these upper limits are higher than the actual planet/star flux ratios. 
Figure [HI shows the upper lir nits together a s a "s pectrum" , and compares it with planet 



atm ospheric models based on iBarman et al.l ( 120051 ) (Barman 2011, private communication) 
and lSudarsky et al.l (120031 ). Our upper limits are at the 6 x 10""^ level on average. The first 
channel gives the best limit and reaches a level of 4.7x10"^, corresponding to a star/planet 
contrast ratio of 2.1x10^:1. This result stands as one of the highest contrast limits achieved 
by closure phase measurements to date. 

Our 90% upper limits for the middle channels are about a factor of 5-8 from the predicted 
value of V And b, suggesting that with further improvement in precision, we will be able 
to start constraining atmospheric models for hot Jupiters with interferometry. In fact, the 
precision of the new calibration for these data is not perfect due to their lack of wide calibrator 
coverages on the sky. In addition, although the new calibration can correct for the closure 
phase drifts caused by dispersion, the use of calibrators from multiple nights makes the night- 
to-night variation hard to calibrate, leaving uncorrected systematic errors in the data. To 
reduce these effects and further improve our calibration precision, the new observing scheme 
using higher spectral resolution (using the grism of R=150 for MIRC in this case) and more 
calibrators is necessary. Benefitting from this experiment and analysis, better precision is 
expected in future observations. 



5. Conclusions and future prospects 

We have simulated the closure phases for the best hot Jupiter candidate v And b, and 
investigated the precision and stability of our measurements obtained with CHARA/MIRC. 
Although our closure phase precisions can reach 0.3°/1.5hrs for the middle wavelength chan- 
nel of MIRC with good conditions, we have also encountered unexpected closure phase drifts 
as large as 8° in other channels. The closure phase drifts are highly correlated with altitude 
and azimuth angles of the targets. The slopes of the drifts also vary across spectral channels, 
centering at the middle channels of MIRC. Because the drifts are correlated for all calibrators 
within an observing run, we are able to model the trend to calibrate the drifts. This new 
calibration model, however, is highly dependent on the altitude and azimuth coverage of the 
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calibrators and may not be accurate for sparse coverage. With a set of diagnostics, we find 
that the closure phase drifts are most likely caused by extra dispersion in the delaylines, 
which contaminates closure phases with non-closed triangles from other wavelengths, espe- 
cially at the edges of the bandpass. Using higher spectral resolution can effectively reduce 
this effect. We therefore advocate future observations of CHARA/MIRC to use the grism 
mode of R=150 for better calibration precision. 

Taking advantage of the well known orbital parameters of v And b, we have combined 
multiple nights of observations together for a joint solution of upper limits for the planet/star 
flux ratios across the H band. Our best upper limit indicates a contrast ratio of about 
2.1xl0'^:l at 90% confidence level, standing as one of the highest upper limits yet achieved 
by closure phase measurements. Future observations with reduced dispersion contaminations 
using the grism mode are expected to have better performance. 

Re cently, photomet ric channels for real-time flux calibration have been implemented for 
MIRC (jChe et al.l |2010| ). The photometric channels have not only improved the visibility 
calibration, but also improved the data taking efficiency for MIRC by a factor of ~ 2 by 
reducing the time spent on reposition ing fibers of each MIRC beam. The CHAMP fringe 
tracker for MIRC (iBerger et al.l |2008[ ) has been commissioned in 2009 and is expected to 
be fully functional soon. CHAMP will help track and stabilize the fringes to increase their 
coherence time, allowing longer integration for MIRC and therefore can increase the SNR to 
the photon-limited regime and signif icantly improve the p recision. In addition, CHARA is 
seeking an adaptive optics upgrade (IRidgway et al.ll2008l ) that could increase the through- 
put of MIRC by at least a factor of 2, which could greatly improve the sensitivity and data 
taking efficiency. With these implemented and upcoming improvements and the new observ- 
ing scheme, we are expecting much better closure phase performance and are optimistic of 
achieving the goal of detecting emission from hot Jupiters in the near future. 
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Fig. 1. — Closure phase simulation of v And b, using the largest triangle of CHARA, Sl- 
El-Wl, and corresponding to an angular resolution of 0.5 mas in the H band. Four out 
of eight wavelength channels of MIRC in the H band are shown. The planet/star flux 
ratios of the four channels (1.52/im, 1. 59Atm, 1.70/im, 1.76/i m) are 0.27x10"^, 0.47x10""^, 
0.6x10"'^, and 0.35x10"'^, respectively (ISudarsky et al.ll2003l ). The orbital parameters used 
in the simulation are listed in Table [2l The four panels assume = 0°, 45°, 90° and 135° 
respectively. The dates of the simulation are chosen arbitrarily. Thanks to the even baseline 
layout of CHARA, the signal level of the closure phase is insensitive to the values of Q, 
although the maximum signal may occur at different times. 
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Fig. 2. — Test observation of v And using the outer array (S1-E1-W1-W2) of CHARA from 
2007Nov22. The top panels show the 60-point averaged data, corresponding to 320 sec of 
integration time. The sohd hnes in the top panels indicate the average closure phase levels, 
while the dashed lines indicate l-cr deviation from the average. The bottom panels compare 
the averaged closure phases with the standard \/N law for normally distributed errors, i.e., 
without systematics. The bottom right panel shows the errors of that channel cannot be 
averaged down beyond 60 averaging points due to the systematic drift. 
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Telescope Triangle: S1-E1-W1 



Slope = -0.059± 0.007 (CP deg/deg) Xv =2.504 
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Fig. 3. — Closure phase vs. Azimuth for six cahbrators in 2008 August. The first, middle and 
the last wavelength channels of MIRC are shown from top to bottom. The data were taken 
with telescope Sl-El-Wl. The red line shows the linear fit of closure phase as a function 
of azimuth. Different colors indicate different targets. The slope of the linear fit and the 
reduced are also labeled in each panel. 
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Fig. 4. — Closure phase vs. Altitude for six calibrators in 2008 August. The red line shows 
the linear fit of closure phase as a function of altitude. Different colors indicate different 
targets. The parameters of the fits are also shown in each panel. The best-fits are good and 
within the scatter and error bars of the data. 
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Fig. 5. — Closure phase as a function of azimuth and altitude using quadratic fit. The 
arrangement and notations are similar to those of Figures [3] and HJ except that the red line is 
a fit to the quadratic plane function: + ai ■ Az + 02 ■ Az'^ + 03 ■ Az ■ Alt + 04 ■ Alt + 05 ■ Alt'^. 
The reduced x^s shown in each panel indicate significant improvement than those of Figures 
Elandll 
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Fig. 6. — The best-fit V^^ of v And witli a uniform disk (UD) model and limb-darkened disk 
(LD) model. The V"^ data are binned into 30 points for better visualization, and are shown 
as filled-dots with error bars. The blue curve shows the best-fit UD model with a reduced-x^ 
of 0.37, obtained from the bootstrap process. The best-fit UD diameter from bootstrap is 
1.097±0.009 mas. The red lii ie slightly below is the LD model with a reduced-x^ of 0.36. 
A power-law limb darkening ( iHestrofferl 119971 ) is assumed for the model. The best-fit LD 
diameter from bootstrap is 1.121±0.007 mas, corresponding to a radius of 1.625±0.011 Rq at 
13.49 pc. 
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Fig. 7. — The closure phases of v And before (left panel) and after (right panel) the new 
calibration scheme. Data from all observing nights are overplotted as a function of closure 
phase vs. wavelength. Due to the variation of wavelength calibration from night to night, 
the data are scattered around each of the eight wavelength channels of MIRC. The dashed 
lines represent the reference level for zero closure phases. The scatters of the data are shown 
in the top right corners. 



-21 - 




0.5 1.0 
Flux Ratio X10"' 





0.5 1.0 1.5 

Flux Ratio x10'^ 



0.0 0.5 1.0 1.5 

Flux Ratio xlO'' 




0.0 0.5 1.0 

Fiux Ratio xlO'^ 



'1.66|xm 




0L£ 

0.0 



0.5 1.0 
Flux Ratio xlO"" 



OL^ 
0.0 



1 .69|xm 



1 .73|xm 



0.5 1.0 
Flux Ratio x10"^ 



OLi. 
0.0 



0.5 1.0 
Flux Ratio xlO"" 



OL 
0.0 



1 .77|xm 



0.5 1.0 
Fiux Ratio x10"^ 



Fig. 8. — Simulated distribution of the best-fit planet/star fiux ratios of v And b in the H 
band, generated from ~150 bootstraps of different nights of data. The dashed fines indicate 
the 90% confidence levels for the upper limits. 
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Barman 201 1 , priv. communication 
Sudarsky et al. 2003 




Fig. 9. — Upper limits with 90% confidence levels for the planet/star flux ratios of v And b 
in the H band, using the newly calibrated CHARA/MIRC data from multiple nights. The 
average upper limit level is 6 x 10"'^. The best channel at 1.52/im shows an uppe r limit 
flux ratio of 4.7x10"^. The solid lines show the latest model based on iBarman et al.l (120051 ) 
(Barman 2011, private communication), assuming a typical radius of 1.3 Rj for the planet. 
The blue line shows the model with incident flux uniformly distributed over the dayside 
of the planet only, while the red line shows the model with full heat redistributi o n ove r 
the entire sphere. The dotted line shows the model prediction from ISudarsky et al.l ( l2003l ). 
assuming a radius of 1 Rj with heat redistribution over the dayside only. 
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Table 1. Observation log 



Observation 


CHARA 


Total integration^ 


Spectral 


Calibrators 


Flux calibration 




Date 


Telescope 


time (min) 


Resolution 


for closure phase 


method 


UT 


2006Oct09 


S2-E2-W1-W2 


8.03 


R=50 


29 Peg, C Per 


chopper 


UT 


2006Octll 


S2-E2-W1-W2 


6.25 


R=50 


29 Peg, C Per 


chopper 


UT 


2006Octl6 


S2-E2-W1-W2 


8.48 


R=50 


29 Peg, C Per 


chopper 


UT 


2007Jul02 


S1-E1-W1-W2 


7.14 


R=50 


7 Lyr, v Peg, a Cyg 


chopper 


UT 


2007Jul04 


S1-E1-W1-W2 


13.38 


R=50 


7 Lyr, v Peg, cr Cyg 


chopper 


UT 


2007Jul08 


S1-E1-W1-W2 


8.92 


R=50 


7 Lyr, v Peg, a Cyg 


chopper 


UT 


2007Jul29 


S2-E2-W1-W2 


3.08 


R=50 


o- Cyg 


chopper 


UT 


2007Jul30 


S2-E2-W1-W2 


3.50 


R=50 


o- Cyg 


chopper 


UT 


2007Aug02 


S2-E2-W1-W2 


10.71 


R=50 


a Cyg, 7 And, 37 And 


chopper 


UT 


2007Aug03 


S2-E2-W1-W2 


9.82 


R=50 


a Cyg, 7 And, 37 And 


chopper 


UT 


2007Aug06 


S2-E2-W1-W2 


26.77 


R=50 


a Cyg, 7 And, 37 And 


chopper 


UT 


2007Augl2 


S1-E1-W1-W2 


40.15 


R=50 


fT Cyg, 7 And, 37 And 


chopper 


UT 


2007Augl3 


S1-E1-W1-W2 


44.17 


R=50 


a Cyg, 7 And, 37 And 


chopper 


UT 


2007Novl4 


S2-E2-W1-W2 


36.58 


R=50 


C Per, a Cyg, v Peg, 7 Tri, 10 Aur 


chopper 


UT 


2007Novl6 


S2-E2-W1-W2 


49.07 


R=50 


C Per, (7 Cyg, v Peg, 7 Tri, 10 Aur 


chopper 


UT 


2007Novl7 


S2-E2-W1-W2 


52.64 


R=50 


C Per, a Cyg, v Peg, 7 Tri, 10 Aur 


chopper 


UT 


2007Novl9 


S1-E1-W1-W2 


29.03 


R=50 


C Per, 10 Aur, 70 Leo, 30 Leo 


chopper 


UT 


2007Nov20 


S1-E1-W1-W2 


30.14 


R=50 


C Per, 10 Aur, 70 Leo, 30 Leo 


chopper 


UT 


2007Nov22 


S1-E1-W1-W2 


52.46 


R=50 


C Per, 10 Aur, 70 Leo, 30 Leo 


chopper 


UT 


20090^22*^ 


S2-E1-W1-W2 


19.63 


R=50 


37 And, 10 Aur, e Cas 


Xchannel"^ 


UT 


2009Oct23'' 


S2-E1-W1-W2 


15.18 


R=50 


37 And, 10 Aur, e Cas 


Xchannel 


UT 


2010Augl3'^ 


S1-E1-W1-W2 


9.30 


R=150 


10 Aur, 7 Tri 


Xchannel 


UT 


2010Augl4<i 


S1-E1-W1-W2 


9.96 


R=150 


37 And, 10 Aur, 7 Tri 


Xchannel 



^Data from each night are split into chunks with less than lOmin of averaging time to avoid smearing of the closure phases. 
''Data taken with linear polarizer 
''Xchannel = photometric channel 

"^Data taJcen with R=150 grism. The new calibration scheme is not applied to these data due to higher spectral resolution. 
For consistency, the 32 spectral channels are averaged into 8 channels in the fits as the other R=50 data. 



Table 2. Adopted Parameters of v And b 



Parameter 


Value 


Ref. 


V (mag) 


4.09 


a 


H (mag) 


2.957 


a 


K (mag) 


2.859 


a 


Distance (pc) 


13.49 


b 


Period (days) 


4.617136 


c 


e 


0.013 


c 


Semimajor axis (AU) 


0.0595 


c 


Semimajor axis (mas) 


4.410 


c 


Tp (JD)A 


2454425.02 


d 


w (deg) 


51° 


d 


Inelination (deg) ^ 


- 58° 


e 


Stellar Diam. (mas) 


1.121 ± 0.007 


f 


Stellar Diam. (R©) 


1.625 ± 0.011 


f 



Re ferences: 
2007^: c. 



a. SIMBAD; b. van Leeuwen 



2006D; 



Wright et all (120091): d iButler et al 



Simpson et al 



(2Ql3); 



f. 

darkened diameter from this work, see 
Figure [HI 



Limb- 
H] and 



^Time of periastron passage 

^Inclination of stellar rotation axis. Note the 
orbital inclination equals to this value only when 
the planet is coplanar with the star. 



Table 3. Results of calibration tests 



Date 


Observing 


Magnitude of slope drift 




Method 


(slopecfts - slopcc/a) 


2008 Aug 22-28 & Sep 02 


Original 


0.37 ± 0.19 


2009 Nov 02 & 04 


Polarizer 


0.13 ± 0.13 


2009 Dec 02, 03, 04 


Polarizer + AQjim slit 


0.18 ± 0.07 


2010 Aug 13 & 14 


Grism R=150 


0.006 ± 0.067 
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Table 4. Calibrator diameters 



Calibrator UD diameter (mas) Reference 



29 Peg 


1.017 ±0.027 


Merand 2008, private communication 


C Per 


0.67 ± 0.03 


getCal^' 


V Peg 


l.Olit 0.04 


Blackwell & Lvnas-Grav f 1994) 


7 Lyr 


0.742 ±0.097 


Leeeett et al. Q9861 


o- Cyg 


0.542 ± 0.021 


Merand 2008, private communication 


7 And 


0.654 ± 0.029 


Based on Kcrvclla & FouQue ('2008) 


37 And 


0.676 ± 0.034 


Based on Kervella & Fouaue C2008") 


7 Tri 


0.522 ± 0.033 


van Belle et al. f2008) 


10 Aur 


0.374 ± 0.079 


van Belle et al. (2008) 


e Leo 


0.677 ± 0.020 


Based on Kervella & Fouquc (2008) 


r] Leo 


0.710 ± 0.033 


Based on Kervella & Fpuquo (2008) 


e Cas 


0.375 ± 0.023 


Based on Barnes et al. fl978) 



''http: / /mscwcb. ipac.caltech.edu/gcWcb/gcWeb.jsp 
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